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Abstract 

We present a Lagrangian-Hamiltonian formalism of a moving dielectric sphere interacting with 
radiation fields. By including the interaction up to the first order in the speed of the sphere, we 
derive the Hamiltonian and perform quantization of both the field and the mechanical motion of 
the sphere. In particular, we show how independent degrees of freedom can be consistently iden- 
tified under the generalized radiation gauge via instantaneous mode projection. Our Hamiltonian 
indicates the form of coupling due to velocity-dependent interactions beyond adiabatic approxima- 
tion. In addition, the Hamiltonian predicts that a geometrical quantum phase can be gained by 
the sphere moving in a light field. 

PACS numbers: 42.50.Wk, 42.50.Tx, 42.50.Pq, 45.20.dc 
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I. INTRODUCTION 



The coupling between an optically levitated sphere and cavity fields has been a subject 



of research interests recently 



-|4(. Such a system corresponds to a basic configuration 
in cavity optomechanics which explores quantum phenomena arising from the interaction 
between mechanical and optical degrees of freedom [5]. In particular, since mechanical 
systems such as a dielectric sphere or moving mirror have masses much greater than that of 
an atom, cavity optomechanics could open a door for testing quantum theory in macroscopic 
scale {4, 6], as well as probing Planck-scale physics Q. Recently, progress in experiments 
las been made in demonstrating cooling and trapping of a dielectric sphere by optical fields 



3] 



In this paper, we present a Lagrangian and Hamiltonian formalism of the interaction of 
a moving dielectric sphere in quantized radiation fields. Previously, Chang et al. [1] and 
Romero-Isart et al. [2] have studied the (adiabatic) interaction Hamiltonian of the system to 
zeroth order in v, where v is the speed of the sphere. Within this accuracy, the interaction 
takes the usual form — | JP ■ E d 3 x (with the polarization P proportional to the electric 
field E) as if the dielectric were stationary. While such a treatment is suitable for a trapped 
sphere at low temperatures, it is interesting to ask for a general theory that takes into 
account interactions due to the motional states of the sphere. In particular, it is known 
that a moving dielectric possesses a velocity dependent polarization and magnetization 
P = (e — 1) (E + v x B) (to first order in v), M = — v x P, and these mot ion- induced 
quantities could modify the dynamics of the sphere and the field. In order to account for 
such effects nonrelativistically, one has to keep the interaction part of the Lagrangian at least 
up to first order in v, and this is our task in this paper - to set up a consistent Hamiltonian 
leading to the quantization of the full system. 

We point out that the moving sphere is a higher dimensional system due to the transla- 
tional and rotational motion in three dimensions, and it requires a treatment different from 

LI 

that in one-dimensional optomechanical systems [9(. For example, an optical field can affect 
the rotational motion of a dielectric via the electromagnetic torque exerting to it 1CH12|. 



and this has been studied in a sequence of experiments [13l4l7l|. Indeed, by including the 
rotation degrees of freedom of the sphere in the Lagrangian, we will show that there is a kind 
of rotational optomechanical coupling caused by the velocity-dependent interactions. There 
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is also some theoretical subtlety in identifying independent degrees of freedom subjected to 
the generalized transverse gauge. Such a gauge is commonly used for field quantization in 



the presence of a stationary dielectric medium [18|-[20|. However, for a moving dielectric, 
the generalized transverse gauge somehow mixes the mechanical degrees of freedom with the 
field. In this paper, we resolve the problem by using instantaneous normal-mode projection, 
and the quantization scheme is compatible with the generalized transverse gauge. 

II. LAGRANGIAN 

The system under investigation is depicted in Fig. 1 in which a rigid dielectric sphere of 
mass m and radius R moves in an electromagnetic field. The sphere is free to rotate about 
any axis through its center of mass (cm.) at q. In our study, we assume a nondispersive and 
nonabsorptive linear dielectric. When the sphere is stationary, its dielectric permittivity is 
given by 

{n 2 , |r — q| < R 
1, otherwise. 

where the convention €q = fiQ = 1 and c = 1 is used for convenience. In addition, the sphere 
is assumed nonmagnetic in its rest frame. 






ielectric sph 



FIG. 1: (Color online) An illustration of an optomechanical system involving a dielectric sphere 
moving freely in electromagnetic fields. The fields can be contained in a generic optical cavity or 
in free space. 

The (nonrelativistic) Lagrangian of the sphere and electromagnetic fields is given by 

L = imq 2 + X -Iu 2 + J d 3 rC(r), (2) 
where / and u> are the moment of inertia and angular velocity of the sphere, respectively. 



The orientation of the sphere is specified by three Euler angles a, (3, and 7 [21], and hence u? 



can be explicitly expressed as u? = (7 sin (3 cos a — (3 sina)e x + (7 sin (3 sin a+$ cosa)e y + (d + 
7COS/3)e z , where e ; (/ = x,y,z) are the basis vectors. £(r) is the Lagrangian density of the 
field after eliminating the electronic degrees of freedom of the dielectric. To obtain the form 
of C, we go to an inertial frame 5"(r) in which the dielectric element at r is instantaneously at 
rest. Assuming the acceleration of the dielectric does not change its macroscopic properties, 
the field Lagrangian density at r in S"(r) is given by the familiar form, 



C' = t ( e E' 2 - B' 2 ) 



(3) 



where E' and B' are the electric and magnetic fields in S'(r), respectively. As the Lagrangian 
density is Lorentz invariant, C can be readily obtained from the Lorentz transformation of 
the fields from S"(r) to the laboratory frame 5*. 

In this paper we confine ourselves to a nonrelativistic motion of the sphere, so that the 
velocity v(r) = q+a; x (r — q) of the dielectric element at any point r satisfies v = |v(r)| <C c. 
By keeping terms up to first order of v in C, the Lagrangian reads 



where 



L = -mq 2 + -Iu 2 + - [ d 3 r (eE 2 - B 2 ) - q • A - u> • T, 
2 2 2 J 

A = J d 3 r(e- 1)(E x B) , 

r = J d 3 r(e - l)(r - q) x (E x B) 



(4) 



(5) 
(6) 



are defined. The Lagrangian 
Salamone 



is a generalization to that of Barton and Eberlein 



22| and 



231 ]. which consider a one-dimensional configuration and focus only on the uniform 



cm. motion of a dielectric slab. Here we will take both q and uj as dynamical degrees of 
freedom which interact with the field through the — q ■ A — uj ■ T term. These velocity 
dependent interaction terms are not considered in previous studies . Interestingly, A 
and r are proportional to the linear momentum and angular momentum of the field inside 
the sphere. 

The Lagrangian should be expressed in terms of the scalar potential V and vector poten- 
tial A via the substitution: E = — dt A — VV and B = V x A. In this paper we shall fix 
the potentials by using the generalized radiation gauge: 



V-[e(r,q)A] = 0. 



(7) 



Such a gauge condition has been employed for performing field quantization in the presence 



of stationary dielectric media 18M20j| . It has the advantage that if the dielectric is at rest, 
the scalar potential V is exactly zero. For a moving dielectric, it can be shown that V under 
condition (J7]) contributes an interaction term of order v 2 in the Lagrangian, and hence the 
effects of V can be consistently neglected in our nonrelativistic Lagrangian which keeps 
interaction terms up to the first order in v. 

It is important to note that the generalized gauge (jTJ) acts as a constraining equation of 
the vector potential A for a given cm. position of the sphere q. Therefore A and q cannot 
be treated as independent degrees of freedom under Eq. (j7j). To overcome the difficulty, we 
identify the independent degrees of freedom via the mode expansion of A: 

A(r,t) = ^Q fc (t)u fc [r,q(t)], (8) 

k 

where {iifc(r, q)} is a set of normal-mode functions of the field when the sphere is sitting 
at rest at the cm. position q. Specifically, the mode function with the mode frequency 
wjfc(q) is defined by 

V x (V x Ufc) - e(r, q)w£(q)u fc = 0, (9) 

subjected to suitable boundary conditions and V ■ [e(r, q)u&(r, q)] = from the gauge 
condition ©. In addition, these modes are orthonormal according to, 

J d 3 r e(r, q)u fe (r, q) ■ Uj(r, q) = 5 kj . (10) 

Note that we have considered as real functions for convenience. For a moving sphere, q 
is a function of time, therefore can be interpreted as an instantaneous normal-mode of 
the field. 

By projecting A onto instantaneous normal-modes, the gauge condition (j7j) is automati- 
cally satisfied and we can treat Qk and q as generalized coordinates, i.e., independent degrees 
of freedom. The Lagrangian (J3J) becomes 

1 1- o 1 



k 



-q-^ikjQkQj - w -^gkjQkQj, (ii) 

kj kj 
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with corrections on the order 0(v 2 ). Here the rj k j and gkj are coefficients depending on q: 
%j(q) = - f d 3 r[e Y ( u k ■ ej) V q (uj ■ e t ) + (e - 1) u fc x (V x uj)), (12) 

l=x,y,z 

gfcj(q) = - J d 3 r{e- 1) (r - q) x [u fc x (V x Uj )] . (13) 

We remark that validity of L in Eqs. (jl]) and (11 ip requires that interaction terms involving 
v 2 and higher orders are negligible. Those v 2 terms omitted in the Lagrangian are corrections 
to translation and rotational kinetic energies such as Smij^qj and SlijUiUj, where Sm^ and 
Slij are respectively the mass tensor and rotational inertial tensor due to the field energy 
stored inside the sphere. Because Sm^ is generally time dependent, it contributes a velocity 
dependent force Sm^qj. Similarly, there is an angular velocity dependent torque SlijCJj 
coming from the rate of change of 81^. If the sphere's motion is consistently described by 
L in Eq. (j3J) to the first order in v, then the force correction SrhijOj and torque correction 
SlijOJj should be small compared with that generated by the — q • A — uj ■ T term. Such a 
condition depends on the physical configuration of the system. In the single mode situations 
with nonzero A and V (see Sec. IV), one can observe that Srhij and Sly oscillate rapidly 
at the field frequency, hence their effects can be averaged out and become negligible in the 
spirit of the rotating wave approximation. 



III. HAMILTON! AN AND QUANTIZATION 

The Hamiltonian defined from L of the system is given by 

H = p-q + HC + ^PkQk-L, (14) 

C=a,/3,7 k 

where p, and tc^ are canonical momenta conjugate to q, Qk, and the Euler angle ( 
respectively. However, since it is rather inconvenient to handle Euler angles, we introduce 
a canonical angular momentum defined by J = (dL/duj), which is related to tt^ via the 
relation: J • u> = air a + fin^ + 77r 7 . From the Lagrangian, p, J, and are given by 

p = mq _ ^2n kj Q k Q j} (15) 

kj 

J = Iu - ^ SkjQkQj, (16) 
kj 

Pk = Qk - fa • *lkj + " ■ Skj) Qy (17) 

3 
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Note that p and J differs from the kinetic momentum mq and angular momentum Iuj for 
nonzero fields, respectively. 

The solutions of q, u> and Qk in terms of the canonical momenta are complicated because 
of the coupled equations ( 1T5|) -( 1T7|) . However, since field-dependent terms that are quadratic 
in v have been discarded the Lagrangian, we have to maintain the consistency by dropping 
terms of same order in eliminating q, uj and Qk for the Hamiltonian. Specifically, we shall 
write 

mq = p + Y^ VkjQkQj ~ P + VkjPkQj, (18) 

kj kj 

Iuj = J + gkjQkQj ~ J + ^2 SkjPkQj- (19) 

kj kj 

By substituting Eqs. ([TBI and ffT9l) into Eq. ffT4|) and again neglecting field corrected inertia 
contributions, we obtain an explicit expression of the Hamiltonian ( TH| . 

k 

where A' = J2kj VkjPkQj and F = J2kj SkjPkQj are defined. 

It is interesting that this Hamiltonian takes a form similar to the minimal-coupling Hamil- 
tonian in electrodynamics, with A' and I" somehow playing the role of vector potential in 
the kinetic energy term. In addition, the Hamiltonian indicates that there is a rotational 
coupling of the sphere described by the second term in ([20]) . It is useful to rewrite T' as 

r' = - J d\ (^) [(r - q) x (H x B)] , (21) 

where B = J2k Qfc(V x Ufc) and II = e J2k Pk u k are the magnetic field and field canonical mo- 
mentum density, respectively. The form of I" is very similar to the field angular momentum 
stored in the dielectric, apart from a proportionality constant. Therefore, approximately 
speaking, the second term of Hamiltonian ([20]) represents an angular momentum coupling, 
i.e., the interaction corresponds to an exchange of angular momenta between the field and 
the sphere. Furthermore, since T'(q) depends on q, there is also a coupling between the 
mechanical rotation and cm. motion of the sphere, mediated by the fields. 

With the classical Hamiltonian ( 1201) . the canonical quantization of the system is readily 
achieved by postulating the dynamical variables into operators with the commutation re- 
lations: [q^pu] = iHS^ [J^J U ] = ihe^ K J K , [Kn,K v ] = -iKe^ K K K , and [Q k ,Pj] = ihS kj , 
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where the Greek subscripts refer to the three axes in rectangular coordinates, and Ku are 



the body-axis components of J 2l|. In this way the quantum Hamiltonian of the system 
takes the same expression as Eq. (I2"U|) . but with P k Qj symmetrized by (P k Qj + QjPk) /2 in 
A' and V. 

In order to represent photon states of the system, we introduce the q-dependent annihi- 
lation and creation operators for each cavity field mode: 



Ofc(q) = A/ 2 ^ fc ( q ) [ukicfiQk + iPh] , (22) 



°* (q)= V2R^) Wq) *- ii5d ' (23) 

which satisfy the commutation relation [a k (q), aj(q)] = 6 k j. Since afc(q) depends on q, for 
each position of the dielectric we have a set of Fock states associated with that position. 
These states can be labeled as | q, where {n k } = {ni,n 2 ,n 3 , . . .} denotes the occu- 
pation number of each photon mode, and £ = (j, m, k) denotes the eigenbasis vectors of J 



(see [2l|): 

J 2 |K},q,0 = ft%' + l)IH},q,£), (24) 

Jz\{n k },ci,0 = frm|{n fc },q,£), (25) 

K z \{n k },q,0 = hk\{n k },q,Z), (26) 

where K z is the z component of J in the body coordinates. Here ({n^}, q, £) is a simultaneous 
eigenstate of the photon-number operator a\.(q)a k (q) and the position operator q i.e., 

4(qH(q)l{ n fc}>q><0 = q>0> (27) 
q|K},q,0 = qlKLq,0- (28) 

Such a set of eigenstates is orthonormal and complete, so that any quantum state of the 
whole system can be expanded in the basis of these eigenstates, i.e., 

I*>=E / rf3 q^(W>q>OIK},q,a (29) 

where C({n k }, q, £) is the probability amplitude. 

With the help of the q-dependent annihilation and creation operators, the Hamiltonian 
Eq. (1201) becomes 

V^'^^ + ^P^^^U^ + l). (30) 



where we have used a shorthand = a&(q) for convenience, and 

ih 



A (q) = - — VkM)]l ~ al a j + % a j ~ a j a k ) , (31) 

. ih \r-^ /Wjfcfq) / + + + + \ 

T (q) = -y ^T^y ( afea i - a k a j + < a i - a >fc ) ■ ( 32 ) 



t„t 



Note that A' and V contains photon-number nonconserving terms a k a, which are responsible 



for photon generation in the dynamical Casimir effect 24j . but this is a subject beyond the 
scope of this paper. For fields at optical frequencies, the aj,a] terms are fast oscillating 
in the interaction picture, and so in the spirit of rotating wave approximation, only the 
photon-number conserving terms a\aj will be kept in A' and V . These terms describe the 
scattering of photons between different modes due to the motion of the sphere. 



IV. SINGLE-MODE SITUATIONS 

In this section, we discuss the Hamiltonian under the single-mode approximation. This 
applies to situations when the field is dominantly contributed by a single mode k, and the 
scattering of photons from the k mode to other modes is negligible within a coherent inter- 
action time. From Eqs. ( 13T|) and ( )32|) . it follows that under the single-mode consideration, 
A' and V only contain photon-number nonconserving terms a£ 2 and a|, and vanish under 
the rotating wave approximation. Hence the corresponding Hamiltonian reads 

p 2 J2 

in which the rotational motion of the sphere is decoupled from the field, and the optome- 
chanical coupling appears only through the position dependent mode frequency w(q). Such 



a form of the Hamiltonian has been considered in Refs. [l, 2]. 

However, we emphasize that the single-mode Hamiltonian (|33l) is based on the real mode 
function in the derivation. The situation is different if complex mode functions are in- 
volved, for example, in a ring cavity which supports traveling wave modes. In the Appendix, 
we show how the Hamiltonian (I5U|) can be modified to incorporate complex modes, in which 
a complex mode function is formed by a linear combination of real modes of the same fre- 
quency. In particular, when photons mainly occupy a complex mode f , by the single mode 
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approximation the Hamiltonian becomes 

,2 



(34) 



where b and w(q) are the annihilation operator and mode frequency associated with the 
f mode, respectively, and A(q) and 7(q) are coupling strengths determined by the mode 
function, 



e £ (r.e,)V q (f.e,) + (e-l)rx(Vxf) 



A(q) = -Mm J d 3 r 
7 (q) = -Mm J d 3 r (e - 1) (r - q) x [f * x (V x f )] . 



l=x,y,z 



(35) 
(36) 



We see that the velocity dependent coupling reappears in the single complex mode. Phys- 
ically, A(q) and 7(q) can roughly be understood as a measure of the field momentum and 
angular momentum stored in the sphere, contributed by a f mode photon. A complex mode, 
such as that of a traveling wave, can carry net field momentum (angular momentum). On 
the other hand, a real mode can be decomposed as a linear combination of complex modes, 
whose contributions of momenta (angular momenta) cancel each other. A familiar example 
is the standing wave mode function sin (k • r), which is a superposition of two traveling wave 
modes e lk r and e~ lk r . This explains why the velocity dependent coupling only appears in 
the complex mode situation. 

To estimate the magnitude of A(q) and 'y(q), we consider a physical situation where a 
dielectric sphere of subwavelength size is placed in a ring cavity, and the field excitation is 
dominantly contributed by a Gaussian beam of an optical tweezer. The configuration of ring 
cavity supports traveling wave modes, so that A(q) and 7(q) can be nonvanishing. Under 
the paraxial approximation, the field mode function is given by, 



f(x,y,z) 



1 Qn 

u(x,y,z)e x + ^7^e 2 



u(x,y,z) = J — 



IT W\Z) 



cxp 



x + y 
w 2 (z) 



exp 



-i tan 1 ( — ) 

zr) 



(37) 
(38) 



where L c is the effective length of the cavity. The beam is linearly polarized in e x , propagat- 
ing along its wavevector k = ke z , with a beam radius w(z) = >j2(z 2 + z R )/kzR (focal plane 
at z — 0), and z R is the Rayleigh range. For a subwavelength sphere satisfying kR <ti 1 and 
R/ zr 1, the Gaussian beam should be a good approximation to the normal-mode, and 
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the q dependence of the mode function should be weak such that the contribution of the 
term involving V q in Eq. f l35]) is negligible, i.e., 

A(q) w -him J d 3 r (e - l)f* x (V x f) . (39) 

Therefore A and 7 can be determined by substituting Eqs. (I3TI) and f[3"5j) into Eqs. fl3"ET) and 
(I3"9"j) . In particular, near the beam focus where their magnitudes are largest, we find 

Ah(n 2 -1) ( R\ 3 2A 

' — (fc^fl) e z , (40) 



3 L c \ZrJ 

4 h(n 2 - l)z R ( R\ 5 n . 2 fi ,ou\flv* q * - \ fA-i\ 
1 = -77 j — (kz R ) (1 + 2kz R ) — y -e x + — e y . (41) 

15 L c \ZrJ \ Zr Zr J 

Numerical calculations were performed based on the parameters in Ref. Q, where n = 1.45, 
R = 100 nm, L c = 4 mm, A = 2n/k = 1064 nm and Zr = 0.53 /xm for the optical tweezer 
at a power P = 15 mW (corresponds to an average photon number (b'b) ~ 10 6 ). We 
find that if the sphere moves under thermal fluctuations at room temperature, the optical 
phase shift U~ l (q • A + uj ■ 7) At accumulated within the coherent time scale At ~ 0.1 ms 
would be on the order 5.1 x 10 -5 rad, which is very small. Furthermore, under the strong 
coherent field of the optical tweezer, we find that the magnitude of the nonadiabatic force 
F = q x (V x A) - V(w • 7) is negligible compared with the restoring force of the optical 
tweezer. 

While the velocity dependent effects on the subwavelength sphere appears to be quite 
weak under the single mode field, we should point out that even under adiabatic motion of 
the sphere, there can be an appreciable effect of geometrical phase due to A. Let us suppose 
that the wave packet of the sphere travels from qj to q^ along a path C, under a strong 
coherent field (e.g., of the optical tweezer). In this case, we may take the field state as 
classical by replacing b^b by its expectation value (b'b). The minimal coupling due to A then 
affects the wave function of the sphere T/>(q) by attaching to it a path-dependent, quantum 
mechanical (geometrical) phase, apart from an overall dynamical phase: 

^( qi ) e-^(q/), Q = h- 1 f dq.A(q)(&t&). (42) 

Jc 



In particular, if the sphere travels along the beam axis of the optical tweezer, 



2, „ ,k 2 R 3 



6 = -r( n ^_i) 



3 V ' L r 



^+(Qz/zr) 2 \Zr 
11 



(43) 



9z 



Hence under the numerical parameters used above, the sphere accumulates a phase of G fa 
— 5.47T by traveling one optical wavelength across the beam focus (i.e. from z = —A/2 to 
z = A/2) (see Fig. [2]). Note that since A is proportional to |f| 2 , the Gouy phase factor (i.e., 
exp[— ita,n~ 1 (z/ zr)}) and the optical phase factor e lkz do not contribute to the calculation 
of 0. 

We point out that in setups which make use of real modes 4], we have 9 = 0, but 
for other configurations involving complex modes, the geometrical phase could be nonzero. 
We also remark that for a subwavelength sphere, G is proportional to the time averaged 
Poynting vector integrated along the path C, which is analogous to the case of an induced 
dipole 25) . However, if the sphere size is appreciable relative to the field wavelength, the 
change of mode structure as the sphere moves should be properly taken into account. In 
our approach, this can be readily achieved by including the contribution of the first term 
involving V q in Eq. (|35|) . 




FIG. 2: Magnitude of geometrical phase |0| accumulated as the dielectric sphere travels from 
q z = —A/2 to q z = A/2 along the beam axis. The beam focus is at q z = 0. We follow the 
numerical parameters as in Ref. [2], where n = 1.45, R = 100 nm, L c = 4 mm, A = 1064 nm, and 
zr = 0.53 fim for the optical tweezer at a power P = 15 mW (hence (b'b) ~ 10 6 ). 
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V. CONCLUSION 



To conclude, we have presented a nonrelativistic Lagrangian-Hamiltonian formalism of 
a moving dielectric sphere interacting with radiation fields. We see that in this three- 
dimensional system, the sphere's cm. degree of freedom q is no longer independent with 
the vector potential A under the generalized radiation gauge, and this poses an interesting 
subtlety in the theory. We have resolved this issue by making use of the instantaneous 
normal-mode projection to consistently identify the independent degrees of freedom subject 
to the gauge, enabling canonical quantization of the system in the usual manner. 

By including the sphere-field interaction up to first order in v, our Hamiltonian (120 j) 
and (I3U|) should capture velocity-dependent optomechanical processes that are not described 
under the adiabatic approximation. For example, coupling between the field (normal) modes 
can result from both the translational and rotational motion of the sphere through the 
interaction described by A' and T', respectively. In addition, these two mechanical degrees of 
freedom become coupled in the presence of radiation fields, due to the coupling characterized 
by r'(q). Such motion- induced coupling can become significant in nonadiabatic regimes, 
especially when the oscillation (or rotation) frequency of the dielectric sphere is close to the 
frequency spacing between two field modes, in which case transitions between the two modes 
can be resonantly enhanced. 

Under the single mode approximation, we have shown that the velocity dependent ef- 
fects are typically very weak for a subwavelength sphere, hence it may justify some of the 
approximations in Refs. Nonetheless, even under such limiting considerations, we 

have also indicated an appreciable geometrical phase (in excess of it) acquired by the sphere 
wavepacket as it moves adiabatically under the single (complex) mode field. Moreover, we 
should emphasize that our theory does not require the single-mode adiabatic approximation. 
With the explicit form of interaction strengths between the sphere and various field modes in 
Eq. ( I30p . our work here enables a further study on the quantum dynamics and applications 
in multimode nonadiabatic regimes. 

Acknowledgment - This work was partially supported by a grant from the Research Grants 
Council of Hong Kong, Special Administrative Region of China (Project No. CUHK401810). 
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Appendix A: The complex mode Hamiltonian 



The complex mode annihilation operators {6fc(q)} are constructed as a linear combination 
of the (degenerate) real mode operators {afc(q)} by 



h{q) = ^2u kj a j {q) 



(Al) 



where the (complex) coefficients Ukj satisfy the unitary property £V Uj-jU^j = J2j UjkU* k , = 
dkk'i so that the commutation relation [&fc(q), &j.,(q)] = 5kw is readily fulfilled. Furthermore, 
since 6fe(q) only mixes degenerate real modes, we have Ukj = except Wfc(q) = u>j(q). The 
complex mode function associated with 6fc(q) is given by 



ffe(r,q) =J3t/^u i (r,q), 



(A2) 



so that in terms of the vector potential A(r) and the field canonical momentum density 

n(r), 

^ 3 r[e(r,q)a; fc (q)A(r) + zn(r)]-f*(r,q). (A3) 



& fc(q) w,, • 

Using Eq. (1A1[) and the properties of U^j, we can rewrite the Hamiltonian ( 13"Uj) as 



= 

where A' and T' reads 
A'(q) = 

r'(q) = 



(p + AQ 

2m 



(j + n 

2/ 



ih 



E 
E 





(q) 




;q) 




(q) 







< ) (q)6 fe 6 j + ^(q^-H.c 



5 g ) (q)6 fc 6 i + gg ) (q)6l6 i -H.c. 



and the coefficients are given by 



^(q) 



»$(q) 



si- } (q) 
gg ] (q) 



/ d 3 r[e ( f * • ej) V q (£, • e { ) + (e - 1) f fc x (V x £,-)], 

l=x,y,z 

J d 3 r[e (4* ■ e,) V q (f, ■ e,) + (e - 1) fj x (V x f,)], 

y d 3 r(e-l)(r-q)x[f fc x(Vxf,)], 
j£ r ( e -l)( r -q)xKx(Vxf,)]. 



(A4) 



(A5) 
(A6) 



(A7) 

(A8) 

(A9) 
(A10) 
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It follows that when the field excitation is dominantly contributed by a single complex mode 



ffc, Eqs. flA5j) and flA6j) can be reduced to, in the spirit of rotating wave approximation 
(where fast oscillating terms such as b\ and b]^ are neglected), 



A'(q) 

r'(q) 



ih r 
T 

ih 



?2(q) 



l8*(q) 



6*6* = A(q)6t&, 



6t6 



(All) 
(A12) 



where the explicit form of A(q) and 7(q) are given in Eqs. (1351) and ( 1361) . and 6 = fefc(q) and 
f = ffc are used for shorter notations. Hence the Hamiltonian (|34|) readily follows in the limit 
of single complex mode. As a remark, a single real mode can (effectively) be decomposed into 
sums of degenerate complex modes whose contributions to the photon-number-conserving 
terms in Eqs. flA5j) and flA6j) cancel out. 
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